A first principles study of the stability and mobility of defects in titanium carbide 



Mikael Rasander,i'2,[*] Biplab Sanyal,i Ulf Jansson, and OUe Eriksson 

^Department of Physics and Astronomy, Division for Materials Theory, 
Uppsala University, PO Box 516, 751 20 Uppsala, Sweden 
^Department of Materials and Nano Physics, KTH Royal Institute of Technology, Electrum 229, 164 40 Kista, Sweden 
'^Department of Materials Chemistry, Uppsala University, PO Box 538, 751 21 Uppsala, Sweden 

(Dated: March 13, 2013) 

We have performed density functional calculations of the formation energies of substitutional 
transition metal (TM) defects, C vacancies, and C interstitial defects in TiC. In addition we have 
evaluated the migration energy barriers for C atoms in the presence of TM impurities. We find that 
the solubility of TM impurities taken from the 3d TM series is low and only Sc and V impurities can 
be dissolved into TiC at equilibrium. In addition, we find that the migration energy barriers of C 
in TiC are greatly affected by the presence of TM impurities: The migration barriers are generally 
lower in the presence of impurities compared to pure TiC and show a clear dependence on the atomic 
size of the TM impurities. We propose that the mobility of C in TiC will be the highest in the 
presence of TM impurities from the middle of the 3d TM series. 



I. INTRODUCTION 

Transition metal carbides (TMC) have for a long time 
been the focus of extensive research. This interest has 
been spawned by their many interesting physical prop- 
erties such as high hardness, high melting temperature 
and high thermal and electric conductivity.^ Properties 
that have given these materials many industrial applica- 
tions, for example as hard coating materials. Recently, 
Wilhelmsson et al.I^lhave shown that by selective alloying 
of TiC with weak carbide forming metals, it is possible 
to grow solutions of ternary TMC that have a built-in 
driving force for the release of C from the carbide system. 
Furthermore, it was proposed that these types of systems 
could be used as low-friction coating materials.^ The al- 
loying reduces the otherwise strong covalent metal-to- 
carbon bonds that are prevailing in TMC, and it is more 
favorable for C atoms to diffuse to the surface by the cre- 
ation of C vacancies in the TMC.'^'^ In experiments on 
thin films of nanocrystalline (Ti,M)C, where the carbide 
phase is embedded in an amorphous C matrix, i.e. (nc- 
(Ti,M)C/a-C), it has been shown that the alloying with 
weak carbide forming metals, M, yields a larger amount 
of C matrix in relation to the carbide phase for systems 
with identical carbon-to-nietal ratio during growth com- 
pared to systems that only contain Ti and C.'^HSl Further- 
more, it has been shown that during heat treatment at 
elevated temperatures there is also an increase of the C 
matrix phase for the alloyed samples which is not present 
in the binary systems.'^"^ When used as a surface coat- 
ing material, the ternary solutions yield a lower friction 
coefficient than the binary sampl es w hile maintaining a 
relatively high degree of hardness.l^El In a previous paper 
we have shown using density functional calculations that 
release of C is a favorable process for solutions of ternary 
(Ti,M)C, where the metal M is any metal from the 3d 
transition metal (TM) series, with the sole exception be- 
ing V where C release never occurs.'^ Furthermore, it was 
found that it is possible to tune the driving force for the 
release of C by careful selection of the alloying metal as 



well as its amount. However, the details concerning dif- 
fusion of C in the ternary solution was not addressed. 

Rock salt TiC is a common and often studied TMC.I^HIll 
Among the 3d TM elements, Ti is involved in the forma- 
tion of the chemically most stable TMC, since the elec- 
tronic structure of TiC is such that the Fermi level is po- 
sitioned in a region with relatively low density of states, 
which separates bonding from anti-bonding states. TiC 
is therefore maximally bonding since the occupation of 
bonding states are complete. By increasing the number 
of electrons in the system, anti-bonding states will begin 
to be filled and the stability of the TMC is lowered. EEMl 
In fact, there are no stable TMC formed by metals found 
late in the 3d TM series. Sc, Ti and V are the only 
3d TM elements that form TMC in the rock salt struc- 
ture, while Cr and Mn only form complicated metal-rich 
TMC similar to Cr23C6.'^Fe, Co and Ni only form meta- 
stable TMC in the cementite structure, i.e. FcsC.^'^' A 
striking feature of TiC is its defect structure since ideal 
stoichiometry is almost never found and t he va riation 
in carbon-to-metal ratio can be rather large.l^^^^^^ The 
mobility of the atoms in TiC has not been extensively 
investigated by theory and especially so for high C va- 
cancy contents. Tsetseris et and Razumovskiy et 
have investigated the migration energy barriers for 
a number of diffusion processes in binary TiC in good 
agreement with experimentally obtained migration en- 
ergy barriers .'^^^'^ll Jt is well established, however, that C 
is the most mobile element in TiC. 

In this paper, we will investigate the stability of TM 
impurities in TiC as well as how the stability of C is af- 
fected by the presence of these impurities. In addition, 
we will address the issue of diffusion of C in TiC in the 
presence of TM impurities. We will show that both the 
stability of C as well as the migration barriers for C dif- 
fusion are greatly affected by the presence of impurities 
in TiC. 

The paper is arranged as follows: In Section [ll] we will 
present the details of the theoretical calculations and in 
Section [Ill| our results will be presented. Finally, in Sec- 
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FIG. 1: (Color online) Crystal structure of rock salt TiC. Ti 
atoms are symbolized by large (red) spheres and the C atoms 
by small (black) spheres. A single 3d TM defect is illustrated 
by the (blue) sphere in the bottom left corner, Ms- Two C 
positions have been labelled by Vi and Vu as is explained in 
the text. In addition, an interstitial C atom is represented in 
the tetrahedrally coordinated interstitial position I. 

tion IIVI we will summarize our results and draw conclu- 
sions. 



II. THEORETICAL DETAILS 

TiC crystallizes in the Bl(rock salt) crystal structure 
which is built up by two intersecting face centered cubic 
(fee) lattices where the Ti and C atoms each occupy one 
of the sublattices. In this study we are mainly interested 
in four kinds of defects which are illustrated in Fig. [T] 
namely (a) single 3d TM defects that substitute for Ti, 
illustrated by in Fig. [l] (b) C vacancies created by 
removing either of the C atoms labelled Vj and Vu in 
Fig. [T] and (c) C interstitial defects illustrated by the 
C atom labelled / in Fig. [l] that occupy the tetrahedral 
position of the metallic lattice. The interstitial defects 
will be set up in two different ways; the first is to insert 
an additional C atom into the tetrahedral site of TiC 
(type I), the second is to move one of the C atoms on 
the C lattice into the tetrahedral interstitial site leaving 
a C vacancy on the lattice (type II). This latter type of 
interstitial defect is more commonly known as a Frenkel 
defect. C vacancies and C interstitial defects have been 
evaluated for pure TiC as well as in the presence of TM 
impurity atoms as is illustrated in Fig. [T] 

The different structures illustrated in Fig. [l] have 
been studied using density functional calculations, where 
we have used the generalized gradient approxima- 
tion (PW91) for the exchange and correlation energy 
functional.'SIl The Kohn-Sham equation has been solved 
using Blochls' projector augmented wave method^^ as it 
is implemen ted in the Vienna ab-initio simulation pack- 
age (VASP).'^^'^ The total energy has been converged 
with respect to the number of k-points, which has been 
set up using the special k-points method of Monkhorst 
and PaclP^ and we have used a plane wave energy cut- 



FIG. 2: Illustration of the diffusion processes that are inves- 
tigated in the present study; Vacancy mediated C diffusion 
(left) and C to interstitial diffusion, i.e. the creation of a type 
II or Frenkel defect (right). The illustrations show a part of 
the (100) surface of the Bl structure with metal atoms rep- 
resented by white and grey circles and C by black circles. 
Black squares represent C vacancies. Grey circles represent 
TM impurities. The arrows show the direction of the motion. 
In case of the C to interstitial process in the right panel the 
interstitial position is represented by the dashed circle. This 
position does not belong to the same plane as the other atoms 
as is clear from Fig. [l] 



off of 600 eV. Relaxations of the atomic coordinates have 
been performed until the Hellmann-Feynman forces act- 
ing between different atoms are smaller than 5 me v/A. 
Spin-polarized calculations have been used when substi- 
tutional defects of Cr, Mn, Fe, Co and Ni have been 
present. However, the magnetic properties will not be 
discussed. 

In order to simulate the mobility of various defects, we 
have calculated the migration energy barriers for diffu- 
sion of atoms in TiC. Here we have considered two types 
of diffusion processes, illustrated in Fig. [2] that are rele- 
vant for C in TiC: The first is vacancy mediated diffusion, 
in which a C atom moves from a position on the lattice 
to a nearby vacant position on the C lattice. The second 
process is related to the formation of the Frenkel type 
of defect, where a C atom moves from the lattice into 
a nearby tetrahedrally coordinated interstitial position. 
The barriers have been calculated using the nudged elas- 
tic band methodj'^SHU! -where we have used a minimum of 
4 images between the different end points. This has been 
tested so that the acquired barriers do not change when 
increasing the number of images. 

The defect formation energies and migration energy 
barriers have been evaluated using the supercell method, 
where the supercells have been built up by 32 primitive 
units cells N, i.e. a total of 64 atomic positions. This 
means that if a single TM impurity atom is substituted 
for Ti, we have an effective concentration of impurities 
on the metal lattice of 1/32 = 0.03125, i.e. about 3 % 
of the Ti atoms are substituted. The same concentra- 
tions are found when evaluating the C vacancy forma- 
tion. These cells are large enough to yield reliable trends 
when it comes to defects formation energies as well as for 
migration energy barriers. Similar supercell size was also 
used by Razumovskiy et al^^ in their study of diffusion 
of intrinsic defects in TiC. 
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Sc Ti V Cr Mn Fe Co Ni Cu Zn 

FIG. 3: (Color online) Evaluated defect formation energies 
for TiC, where Esoi, Eg, and By have been evaluated accord- 
ing to Eqs. ([2|, ([3|, and Q respectively. The labels on the 
X-axis marks the TM impurity. Note that Ti represent the 
case of pure TiC. 

III. RESULTS 

The ability for TM impurities to be dissolved into TiC 
is guided by the reaction 

(7V-l)TiC + M, ^TiA,_iM,C^_i, (1) 

where A'^ — 1 is the amount of TiC and Mg is an impurity 
atom. The expression on the right of Eq. ([T]) is equivalent 
to having a TM atom Ms dissolved into TiC by the cre- 
ation of Tiftr-iMsCAT-i, i.e. a TM atom is positioned on 
the Ti lattice followed by the expansion of the lattice by 
the formation of a C vacancy. We therefore evaluate the 
energy for dissolving TM impurities into TiC according 
to 

Esol - £;(Tiw-iM,CAr_i) - ^(TiAr_iCAr_i) - S(M,), 

(2) 

where i?(Tijv_iMsCAr_i) is the total energy of TiC with 
one substitutional 3d TM defect M^ and one C vacancy, 
-E(Tijv-iCAr_i) is the energy of defect free TiC, and 
E{Ms) is the energy of the TM in its standard refer- 
ence state.'^ The result is shown in Fig. [3j where it is 
found that Sc and V are the only impurities that may be 
favorably dissolved into TiC. This is because Sc and V 
arc the only metals in addition to Ti that form a TMC 
in the rock salt structure and our result is in agreement 
with the results obtained for alloyed solutions of (Ti,M)C 
in Ref. 3, 

We also define the energy required to substitute a TM 
atom for Ti by 

Es = E{TiN-iM,CN) - E(TiNCN) - {E{Ms) - E{Ti)} , 

(3) 

where £'(TiAr_iMsCiv) is the total energy of TiC with 
one substitutional 3d TM defect M^, E{TiNCN) is the 
total energy of defect free TiC, and E{Ms/Ti) is the en- 
ergy of the metal M^/Ti. In addition, the formation en- 
ergy of a C vacancy is given by 

E^ = E{TiN-iMsCN-i)-E(TiN-iM,CN)+E{C), (4) 



where the notation is consistent with Eqs. ([2]) and 
and E{C) is the energy of one C atom in a graphite unit 
cell consisting of 4 C atoms. By the use of Eqs. ^ and 
Q it is possible to rewrite Eq. ^ as 

Esoi=Es+E^ +E^ (TiC), (5) 

where £^-'^(TiC) is the formation energy of TiC, i.e. 
E^{TiC) = E{TiC) - E{Ti) - E{C). As is evident from 
Eq. ([5]), the ability to dissolve TM impurities into TiC 
depends on the energy required to substitute a Ti atom 
in TiC for a TM impurity (Eg), the ability for C vacancy 
formation in the presence of a TM impurity (E^) as well 
as the formation energy of TiC. The energy required to 
substitute a TM for Ti and the formation energy of a C 
vacancy are also shown in Fig. |3] It is clear that the sub- 
stitution of a TM for Ti is an unfavorable process with a 
very high formation energy, especially towards the end of 
the TM series. This behavior is understood by consider- 
ing a rigid band model where the substitution of Ti for a 
TM atom with a larger number of valence electrons leads 
to the occupation of antibonding states which increases 
the total energy of the system compared to pure TiC, 
see the discussion in Section [l] In the case of a Sc impu- 
rity, the occupation of binding states is lower compared 
to the pure TiC which also increases the total energy of 
the system. The C vacancy formation in the presence of 
TM impurities is only unfavorable in the case of Sc and V 
impurities. A transition is made for a Cr impurity after 
which the C vacancy formation energy becomes nega- 
tive. This means that in the presence of these latter TM 
defects, C vacancies will spontaneously be formed. In 
addition, we note that in the presence of any of the TM 
defects, the C vacancy formation is lower than in pure 
TiC, which means that the binding strength between the 
TM and the surrounding C atoms is weaker than between 
Ti and C. 

We conclude that the main reason for the difficulty of 
dissolving TM impurities into TiC stems from the substi- 
tutional energy (Eg) component of Eq. ([s]). The energy 
required for the substitution of Sc and V for Ti is low 
enough to be overcome by the formation energy of TiC, 
resulting in a negative solubility energy Esoi- For the 
other TM defects, the substitutional energy is so large 
that the strong binding in TiC is not enough to allow 
for dissolving these metals into the system, even though 
the C vacancy formation energy becomes negative for the 
later TM elements. 

In addition, we note that the C vacancy formation en- 
ergy, given by Eq. Q, corresponds to the formation of 
a C vacancy by the removal of a C atom in TiC that 
is then moved out of the carbide and into the form of 
graphite. This reaction is valid when there is graphite 
present with which C can be exchanged, e.g. for the 
case of a carbide phase embedded in a C matrix. Such 
a system is considered to be rich of C. For Ti rich (C 
poor) conditions, C vacancies can be formed when Ti is 
dissolved into the TiC phase, described by the reaction 
in Eq. ([I]) with Mg = Ti and evaluated using Eq. ([2|. 
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FIG. 4: (Color online) The calculated C vacancy mediated 
migration energy barriers Em in the presence of TM impuri- 
ties. Note that Ti represent the case of pure TiC. The inset 
shows the metallic radius R of the TM atoms. The metallic 
radii given here are taken from Ref. 1371 



Ur 



12- 



10- 



® 6 - 



2- 



□ e; 
■ e; 

E., 




Sc Ti V Cr Mn Fe Co Ni Cu Zn 



FIG. 5: (Color online) The calculated formation energies of 
type I and II C interstitial defects according to Eqs. ([6| and 
TI). The migration energy required to from the type II defect, 



-Ejvl is also shown. Note that Ti represent the case of pure 
TiC and, in addition, the migration barrier in the presence of 
Sc impurities is not evaluated. 



This energy is shown in Fig. |3j and it is negative which 
means that for Ti rich conditions, Ti will be dissolved 
into TiC with the associated formation of C vacancies. 
This shows that it is more favorable to have C vacan- 
cies in TiC than it is to form a separate Ti phase and 
can be related to the fact that sub-stocihiometric TiC^: 
with .5 < a: < 1 is stable against decomposition into Ti 
and C.l^*^ We note that the difference between the C va- 
cancy formation energy for Ti rich conditions, evaluated 
according to Eq. ^ , and the C vacancy formation energy 
for C rich conditions, evaluated by Eq. ([4|, is identical to 
the formation energy of TiC, E^(TiC), as can be seen in 
Eq. (|5|, since in the case of pure TiC the substitutional 
defect formation energy Es = 0. C vacancy formation 
in pure TiC for Ti rich conditions is therefore favorable, 
while it is unfavorable for C rich conditions. Even so, the 
evaluated trends in the vacancy formation energy in the 
presence of TM impurities will vary in the same fashion, 
irrespective of the external conditions. 

In Fig. |4] we present the evaluated vacancy mediated 
migration energy barriers of C in the presence of TM im- 
purities. The value obtained for this process in pure TiC 
(second staple from the left in Fig.|4| is in agreement with 
the value obtained by Tsetseris et of 3.46 eV. The 
difference between our value and the value of Tsetseris 
et al. is attributed to the use of different sized super- 
cells. Our study is done using a smaller supercell, which 
will yield a larger value for the migration energy barrier 
compared to a larger supercell, since atoms are allowed 
to relax further in a larger supercell and thereby lowering 
the energies along the migration path. The migration en- 
ergy barrier also correspond rather well with experimen- 
tal observations (3.30 eV,E2l 3.60 eV, 1231 4.15 eWH). We 
note, however, that lower barriers have been obtained ex- 
perimentally in the order of 2.44-2.70 eVl^^lbut these have 
been related to a different diffusion process by Tsetseris 
et al.li^ where the C atoms migrate between interstitial 
sites and not between lattice sites on the C lattice as has 



been evaluated here. 

According to the data shown in Fig. |4] the vacancy 
mediated process has a higher barrier if it takes place in 
the presence of Sc impurities compared to pure TiC. All 
other TM impurities yield a smaller barrier. We also note 
a clear trend along the TM series: The highest values are 
found in the beginning and the end of the series, with a 
clear minimum in the migration energy barrier in the 
middle. This trend does not reflect the same behavior as 
the stability of the TM impurities, neither does it reflect 
the stability of the C atoms in the vicinity of the TM 
impurities, see Fig. [3] Rather, we find that the behavior 
of the migration energy barriers reflects the behavior of 
the size of the TM involved. In the inset of Fig. |4j we 
show the metallic radius of the TM atoms. As can be seen 
in Fig. |4j the metallic radius along the 3d TM series has 
its largest value for Sc, with a minimum in the middle of 
the series, and a low rise at the end. The finding that the 
height of the migration energy barriers is related to the 
atomic sizes is reasonable since a large impurity atom will 
inhibit the motion more than a smaller impurity, which 
is exactly what is shown in Fig.|4j In the past, the size of 
the TM impurities have also been found to correlate with 
other materials properties, e.g. the tendency for surface 
segregation of TM impurities of the TiC(lOO) surface.l^ 

We now turn the attention to C interstitial defects and 
the results are shown in Fig. |5] The interstitial defect of 
type I in the presence of TM impurities is evaluated as 

£;/ = E{Tin-iM,Cn+i) - (Tijv-iM.CAr) - £;(C), (6) 

where i?(Tijv_iMsCAr+i) is the energy of a system with 
an additional C atom in the interstitial position / in 
Fig. [T] while the interstitial defect of type II is evalu- 
ated as 

£f = E{TiN-iM,C*j^) - E(TiN-iM,CN), (7) 

where i?(Ti7v-iMsC^) is the energy of a system with a 
C atom at the interstitial position / and a vacancy at 
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position Vi in Fig. [T] We find that the formation energy 
of a type I C interstitial defect (shown in Fig. [5]) has a 
minimum in the presence of V and Cr defects and that 
this energy becomes very large towards the end of the 
3d TM series. This behavior is also related to the size 
of the TM elements, but also strongly related to the TM 
elements ability to form TMC which is strongly reduced 
towards the end of the TM series. The type I defect has 
an additional C atom in the tetrahedrally coordinated 
interstitial site next to a TM atom and it is not surprising 
that the formation energy is high and that it shows a 
clear atomic size contribution since the density of atoms 
is large. The type II defect, or Frenkel defect, on the 
other hand has a much lower formation energy. In Fig. [5] 
we see that the type II defect is easier to form close to 
TM impurities compared to the pure TiC and towards 
the end of the TM series the formation energy becomes 
rather small. 

In Fig. [5] we also show the migration energy barriers 
for moving a C atom on the lattice, Vj in Fig. [l] to the 
interstitial site / in Fig. [l] We find that this barrier has a 
similar behavior in the presence of TM impurities as the 
C vacancy mediated migration barrier shown in Fig. [4j 
which is again related to the sizes of the TM impurities 
compared to Ti. 

IV. SUMMARY AND CONCLUSIONS 

In summary, we have performed density functional cal- 
culations of the formation energies of TM impurities, C 
vacancies, and C interstitial defects in TiC. We find that 
the substitution of TM elements for Ti is an unfavorable 
process and, in fact, the only elements that can be fa- 
vorably dissolved into TiC at equilibrium are Sc and V. 



However, it has been found that during non-equilibrium 
conditions large quantities of TM, e.g. Fe, Ni, Ft and 
Al, can be dissolved into TiC to form ternary TMC 
systems.l^^l^ Furthermore, we find that the presence of 
TM impurities greatly affects the stability and mobility 
of C atoms in the vicinity of the TM. Especially, the C 
vacancy formation energy is always lower in the presence 
of the impurities compared to pure TiC, and therefore 
the impurities will attract C vacancies. We also find that 
the migration energy barriers of C are always lower in the 
presence of TM impurities compared to pure TiC, except 
for the case of Sc, due to the atomic size effect of the TM 
atoms compared to Ti: If the size of the TM is smaller 
than Ti, the migration energy barrier will be lower. This 
means that the mobility of C is higher in systems with 
large quantities of TM impurities, and especially so for 
TM impurities from the middle of the 3d TM series. We 
therefore propose that the C release effect shown to exist 
for alloyed solutions of (Ti,M)Cl2HIIwill be more favorable 
if the alloying element is taken from the middle of the 3d 
TM series. 
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